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An overview is given of the limitations of Luttinger liquid theory in describing the real 
time equilibrium dynamics of critical one-dimensional systems with nonlinear dispersion 
relation. After exposing the singularities of perturbation theory in band curvature effects 
that break the Lorentz invariance of the Tomonaga-Luttinger model, the origin of high 
frequency oscillations in the long time behaviour of correlation functions is discussed. 
The notion that correlations decay exponentially at finite temperature is challenged by 
the effects of diffusion in the density-density correlation due to umklapp scattering in 
lattice models. 



1. Introduction 

When Haldane coined the name "Luttinger liquid" JD the point was to emphasize the 
universality of the theory beyond the exactly solvable model studied by Tomonaga 
and Luttinger PES Indeed, the mapping from interacting fermions to noninter acting 
bosons that renders the Tomonaga-Luttinger (TL) model solvable depends crucially 
on the approximation of a linear dispersion relation for low energy excitations. 
However, the thermodynamic properties predicted by Luttinger liquid (LL) theory 
are asymptotically exact in the low energy limit for generic critical one-dimensional 
systems because perturbations associated with band curvature are irrelevant in 
the renormalization group sense. This explains why Luttinger liquid behaviour is 
observed in so many different systems, such as quantum wires, carbon nanotubes, 
spin chain compounds and cold atoms in optical lattices.^ 

Neglecting irrelevant perturbations, the elementary excitation of the one- 
component TL model is a free boson with massless relativistic dispersion u q = v\q\. 
As a result, the model is Lorentz invariant: Its correlation functions are preserved 
by continuous rotations between space and Euclidean time. More precisely, the TL 
model is invariant under local conformal transformations in (1 + 1) dimensions. 
This allows one to bring in the arsenal of conformal field theory (CFT), with cen- 
tral charge c = 1, to compute correlation functionsP Conformal invariance dictates 
that for large distances x and long real times t the ground state correlation function 
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(*(x,t)*H0,0)) ~ - 7 -— -— - , (1) 



for a given field decays as a power law 

1 

(x - vt) 2A ~x + vtf 

where A + and A_ are conformal dimensions. The conformal dimensions of physical 
operators are determined by the Luttinger parameter K of the TL model. Both 
the velocity v and the Luttinger parameter K can be extracted from the finite size 
spectrum in general, or calculated exactly for integrable models in particular P The 
result can be generalized to models with more than one gapless degree of freedom, 
such as the Hubbard model away from half-filling P Using the conformal mapping to 
the cylinder geometry with compactificd Euclidean time direction, one can predict 
that at finite temperatures correlation functions decay exponentially in both x and 
tJffl " 

Perhaps due to the remarkable success of CFT methods in LL physics, the valid- 
ity of Eq. (JTJ is often overstated. The truth is, for any model in the LL universality 
class where the dispersion relation is not exactly linear, the CFT result for time- 
dependent correlation functions does not give the correct long time behaviour for 
real time t > \x\/v, i.e. inside the light cone. While there are examples of gap- 
less one-dimensional systems to which CFT techniques clearly do not apply — for 
instance spinful fermions in the spin- incoherent Luttinger liquid regime^ or the fer- 
romagnetic Bose gas^l — the general reason for the breakdown of LL theory in 
real time dynamics is the effect of band curvature. Although formally irrelevant, 
perturbations to the TL model that take the form of boson decay processes gener- 
ate singular contributions to dynamical correlation functions. The solution to this 
quandary gave birth to the subject of "nonlinear Luttinger liquids"^ (see Ref. 12 
for a detailed review). Some predictions of the field theory for nonlinear LLs have 
been recently confirmed by an exact form factor approach P^l 

The purpose of this chapter is to provide an overview of the breakdown of LL 
theory in equilibrium dynamics and discuss its consequences for the long time be- 
haviour of correlation functions in critical one-dimensional systems. We will mainly 
focus on two aspects: the contribution of high energy modes to correlation functions 
at zero temperature and the diffusive contribution due to umklapp processes which 
dominates the long time tail at finite temperatures. 



2. Breakdown of Luttinger liquid theory by band curvature effects 

Fermi liquid theory breaks down in one dimension because scattering between 
two disconnected Fermi points always leads to singularities in the particle-hole 
and particle-particle channels, making quasi-particles unstable as the fermion self- 
energy diverges.^ Similarly, LL theory fails to describe dynamical response func- 
tions because the boson self-energy due to band curvature terms is singular. How- 
ever, the singularity in this case is connected with the macroscopic degeneracy of 
states comprised of multiple bosons with the same chirality, which is an artifact of 
the linear dispersion approximation in the TL model. 



July 23, 2012 0:58 WSPC/INSTRUCTION FILE LLpereira 



Instructions for Typing Manuscripts (Paper's Title) 3 

Consider the TL model for spinless fermions (following the standard 3-ology 
notation^ 



-R,L 



(2) 

Here r = R, L = ± denotes right and left movers, defined from single-particle states 
with momentum around ifcp, and :: refers to normal ordering with respect to the 
free fermion ground state. The fermion field operator reads &(x) e fx iPr(x) + 
e~ fX iI)l(x). Bosonization 4 maps the fermionic model to the Gaussian model 

H = E / dx \ [ ; ( 5 -^) 2 : + : ( a ^) 2 ; ] > ( 3 ) 

r=R,L t ' 

where the normal ordering is with respect to the vacuum of bosons and the chiral 
bosonic fields obey the commutation relation [ip. r {x), d x /ip r /(x')] = —ird r y5(x — 
x'). For (72,34 <C Vp, the renormalized velocity is given approximately by 
v « vf — 34/(27r) + 0(g 2 ). In the notation used here, the chiral fermion fields are 
bosonized as 

tp r (x) « (27ra)~ 1/2 exp[— i\/2n(X(p r + \<P- r )], (4) 

with A = (y/K + l/\/K)/2 and A = (-s/K + l/VK)/2. Here a is a short- 
distance cutoff and K is the Luttinger parameter given approximately by K &s 
1 — 92/(2ttvf) + 0(g 2 ). Within the TL model, the single fermion Green's function 
at zero temperature is given by G(x,t) = -i(^(x, t)* t (0, 0)) = e lkFX G R (x,t) + 
e~ lkFX GL(x,t), where the Green's functions for the chiral fermions reacP^ 

G r {x,t) = -i(i/v(a:,i)V4( : )) = (r/2n)(x - rvt)- x2 {x + rvt)~^ . (5) 
The result in Eq. §5§ has the form of Eq. (1) with conformal dimensions A± = 
| + 4± 2r) . Taking the Fourier transform with the proper time ordering pre- 
scription one finds that the particle addition part of the single- fermion spectral 
function does not have a quasi-particle peak; instead it behaves as a power law 
above a threshold energy, A{k, u) ~ k x (w — v5k)~ 2+x for 8k = k — kp -C kpW&^ 
Importantly, the exponent in the spectral function is directly related to the singu- 
larities of G r (x,t) in Eq. ([5]) along the light cone x — ±vt. For 32,34 "C vp, the 
interaction-dependent exponent can also be obtained by resumming logarithmic di- 
vergences in the perturbation theory in the fermion description,^ but the advantage 
of bosonization is that the fermion interactions are treated easily by the rescaling 
of the bosonic fields leading to Eq. (|4j . 

What if we perturb the TL model with band curvature effects, which are indeed 
present in any real system? Let us assume, as usual, that we are allowed to truncate 
the Hilbert space to low energy states around the Fermi points, but now we add a 
parabolic term in the dispersion^ 

5H = --!- f dx : ( 6 ) 



-R.L 
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where m is the effective mass at the Fermi level. After bosonizing and performing 
the Bogoliubov transformation that diagonalizes the TL model, the band curvature 
term generates two types of operators in general 



SH=^ Jdx {ry_[: (d x <p L ) 3 : - : (d x cp R ) 3 :] 

+77+ [: {d x ip L ) 2 d x ipR : - : (d x ip R ) 2 d x <p L ) :]} , (7) 

where rj± are coupling constants of order 1/m which can be calculated to lowest 
order by bosonization or fixed by phenomenological relations! 19 ^ 20 ! The important 
point is that these cubic terms spoil the solvability of the TL model since they 
introduce interactions between the bosonic modes. Nevertheless, we might hope they 
can be treated perturbatively. To illustrate the problem with perturbation theory, 
it suffices to consider the ?y_ term, which does not mix the two chiral components of 
the bosonic field. Calculating the first-order correction to the single-fermion Green's 
function, one finds^ 



5G r (x,t) iA 3 ??- 
G r (x,t) v 



1 x + rvt 



1 x — rvt 



+ rvt (x + rvt) 2 



x — rvt [x — rvt) 2 

(8) 

For x — > oo or t — » oo, the expression in Eq. ([5]) decays faster than the result for 
the TL modelEI This is expected since simple power counting tells us that the rj- 
perturbation has scaling dimension three and is irrelevant at the LL fixed point. 
Hence the argument for the universality of LL theoryP If we are interested in low- 
energy thermodynamic properties, band curvature effects are harmless because they 
only give subleading corrections to correlation functions at large distances x 3> a, vt. 

The problem shows up when we are interested in dynamical response functions 
which depend on both momentum and frequency, such as the spectral function 
A(k,u>). Remarkably, the correction in Eq. (|8|) is more singular at the light cone 
x = ±vt than the unperturbed LL result in Eq. ([5]). For the spectral function, this 
means that the corresponding correction 5A(k,uj)/A(k,uj) ~ r]_5k 2 /(uj — vSk) has 
a singular frequency dependence and diverges more strongly at the lower threshold 
uj = v5k than the LL result. The singularity actually gets worse at higher orders of 
perturbation theory in ?7_, as one obtains more powers of l/(x ± vt) in 5G r (x, t). 

The same problem is present in the calculation of the density-density correlation 
function x(x, t) = {n(x, t)n(0, 0)). In the boson description, the fluctuation of the 
fcrmion density operator is represented by 

n(x) = : * t (a;)*(a;) : « y/Kj%d x <t>{x) - (l/2na) cos[V '4ttK<P(x) - 2k F x], (9) 

where <f>(x) = [ipl(x) — ipr(x)]/v2. The correlation function for the long-wavelength 



a Actually, the correction due to r\— vanishes for t = 0. Equivalently, the correction for equal-time 
correlation functions calculated around Eq. (5.4) of Ref. 1 vanishes if the "angle" of the Bogoliubov 
transformation that mixes right and left movers is set to zero. 
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part of n(x) is equivalent to the boson propagator 

X (x,t) ~ V<: ^Vv : (x,t) : : (0,0)) = -<^(x,t)9 x 0(O,O)). (10) 

* — ' 7T 

r,r' 

The Fourier transform of x( a; : f° r small momentum q <^ kp yields the dynamical 
structure factor S(q, iS) = — 2lmxret(q, w) ps If v\q\). The delta function peak 

obtained as an approximation for S(q -C kp,uj) within the TL model corresponds 
to the spectral function of the coherent bosonic mode with well defined energy and 
momentum. Following the analogy with Landau's Fermi liquid theory, one would 
expect that going beyond the TL model and introducing band curvature effects 
would lead to a finite boson lifetime. However, the attempt to calculate a boson self- 
energy using perturbation theory fails J2Sl Similar to the fermion Green's function, 
perturbation theory in rj- for x( x i generates corrections which are increasingly 
more singular at the light cone x = ±vt. To understand this singularity, we note 
that in momentum and frequency domain the three-legged vertex ?7_ in Eq. ([7]) 
allows the single boson with momentum q to decay into two bosons with momenta 
qi and <?2 = q — qi (see Fig. BJa)); however, the energy of the intermediate state 
when the boson lines are put on shell is always w = vqi + vqi = vq, independent 
of the internal momenta. This huge degeneracy is present at any finite order of 
perturbation theory. The result is that to any finite order the boson decay rate due 
to the r\- perturbation diverges on shell as ~ 6(oj — vq) . 




Fig. 1. (a) In the boson description, band curvature is represented by an interaction vertex 
that scales with the momenta of the bosons. But due to Lorentz invariance of the unperturbed 
TL model, the two-boson state after the decay is degenerate with the initial single-boson state, 
regardless of the value of qi (or 52 = 9 — 91 )• (b) The divergences of perturbation theory in 
band curvature operators can be resummed by refermionizing to new fermions with renormalized 
nonlinear dispersion. Density excitations are now represented by non degenerate particle-hole pairs, 
(c) Broadened boson peak in the density structure factor S(q, u)) of a nonlinear Luttinger liquid 
for q <^kp. 
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The degeneracy of the many-boson states stems from the linear dispersion ap- 
proximation of the TL model. In the model of free fermions with nonlinear disper- 
sion, particle-hole pairs with the same total momentum — a linear combination 
of which defines the bosonic excitations — are not degenerate. Thus the cause of 
the breakdown of LL theory is the impossibility of starting from the TL model and 
breaking Lorentz invariance within finite order perturbation theory. The picture of 
a "quasi-boson" with well behaved self-energy due to band curvature simply does 
not work. It is frustrating that the bosonization method, so helpful in resumming 
the divergences of Fermi liquid theory in one dimension, cannot handle the innocent 
looking perturbation in Eq. ([6]), which is quadratic in fermions. 

The question then is whether there is an alternative representation which cap- 
tures the essential physics and resums the divergences of band curvature effects in 
LL theory. Fortunately, the answer is yes. The trick is to refermionize the Hamilto- 
nian with the cubic perturbation rj- to obtain a model of noninteracting fermions 
with nonlinear dispersiorP^H 

H + SH^ ]T fdx; 4>l [v(-ird x ) + ^t-(-td x ) 2 ] j r : (11) 

r=R,L~' 

The new fermions are defined such that : V4Vv := —rd x (p r /\/2ir and differ from the 
original fermions by string operators We note that the fermion interactions (72,54 
in Eq. ^ are absorbed into the renormalization of the velocity v and effective mass 
ry_ . Including the ry_ operator as the parabolic term in the dispersion of the new 
fermions lifts the degeneracy of the many-boson intermediate states. The approxi- 
mation in Eq. ((TTj) is to neglect the 77+ operator defined in Eq. ([7]) as well as more 
irrelevant operators (with dimension four and higher). Within this approximation, 
the free-boson peak in S(q, u) broadens into a two-fermion continuum with rectan- 
gular line shape and width 5oj q ~ ^-q^P^IlThis result is exactly what one expects 
from summing the infinite series of diagrams in the r;_ perturbation.^ 

The renormalized fermion band with nonlinear dispersion relation offers a conve- 
nient starting point to study dynamical response functions. With Lorentz invariance 
broken at the outset, the thresholds of the exact spectrum are not the same as in LL 
theory. While the rj- band curvature term sets the width of the two-fermion con- 
tinuum for q <^ kF, additional irrelevant operators generate interactions between 
the new fermions and give rise to power law singularities at the edges of the spec- 
trum. Rather than governed by light cone effects, the new singularities of nonlinear 
Luttinger liquids are in analogy to the x-ray edge problem of optical absorption 
in metalsP^l The lower threshold u;_ (q) below which S(q, w) vanishes is given by 
the minimum energy of a particle-hole excitation with momentum q created in the 
renormalized fermion band. For positive band curvature ry_ > as in Fig. 1(b), this 
corresponds to a particle at the Fermi point and a hole as deep as possible with 
energy w_ (q) — vq — r/^q 2 /2. Like in the x-ray edge problem, this "deep hole" can be 
described as an effective quantum impurity which propagates with different velocity 
than the low energy modes. The power law singularity at the edge ^J-(q) has an 
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exponent proportional to q in the case of short range interactions.^ Moreover, in- 
teractions between right and left movers make the spectral weight extend above the 
upper threshold of the two-fermion continuum uj + (q) = vq + r]-q 2 /2 predicted by 
the approximation in Eq. pip and S(q, uj) acquires a tail that decays as ~ rf^q /u> 2 
at high frequencies w — vq ^> r7_ g 2 E2H25] 

The picture that emerges for the broadening of the peak in S(q <C kp,uj) due to 
irrelevant operators is illustrated in Fig. 1(c). Unlike the Lorentzian quasiparticle 
peak in the spectral function of Fermi liquids, the "quasi-boson" peak of nonlinear 
Luttinger liquids is asymmetric and has an x-ray edge type singularity above the 
lower threshold w_ (g) . 

Accounting for band curvature effects, the singularities of the single fermion 
spectral function of Luttinger liquids are also modified!^ In LL theory particle and 
hole Green's functions coincide because the TL model is particle-hole symmetric, 
but in the presence of band curvature this is no longer the case. For positive band 
curvature, the excitation that creates a single deep hole defines the lower thresh- 
old of the support of the hole spectral function for a given momentum fc. As a 
result of kinematics, the power law singularity at the deep hole threshold cannot 
be broadened by any interactions. In contrast, the support of the particle spectral 
function does extend below the energy of the single particle excitation. In this case, 
three-body scattering processes in generic (i.e. nonintegrable) models allow the sin- 
gle particle to decay into the continuum and a Lorentzian peak with decay rate 
7^ oc (fc — kp) 8 is obtained.^ The only surviving power law in the particle spectral 
function has a positive exponent and is found at the absolute lower threshold of the 
support, located at the energy of the deep hole excitation. 

Going back to real space and time, we should expect the long time decay of 
correlation functions of Luttinger liquids to be strongly affected by the x-ray edge 
type singularities at the thresholds of the nonlinear spectrum. This will be the 
subject of the next section. 

3. High energy contributions to time-dependent correlation 
functions 

Once Lorentz invariance is broken, the exponents for the long time decay of corre- 
lation functions are not constrained to be the same as the ones for large distance 
decay. In fact, contrary to conventional wisdom, the long time behaviour at zero 
temperature is not even dominated by low energy modes, but by high energy saddle 
point contributions which take advantage of the dispersion nonlinearity. 

To see how this comes about, consider the simple case of free fermions with 
parabolic dispersion relation e(k) — k 2 /2m. The hole Green's function is given 
exactly by 

G h ( X ,t) = (*t(a;,i)*(0,0)) = / " — e ^^(k 2 -kl)t/2rn^ ( 12 j 

J-k F 27r 

The real part of Eq. (fT2|) is plotted in Fig. 2(a). There are clearly two distinct regions 
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in the (x,i) plane. Outside the light cone, |x| > vpt, the Green's function oscillates 
with distance, but not with time. This is consistent with the usual contribution 
from the low energy modes which come with factors of e ±lkFX . For large distances, 
\x\ 3> hp , vpt, the power law decay of Gh(x,t) is well described by Eq. flSJ with 
K = 1. In contrast, outside the light cone, \x\ < vpt, the Green's function exhibits 
time oscillations in addition to spatial oscillations. Such time oscillations can only 
come from modes with finite frequency. Indeed, a moment's reflection shows that 
the integral in Eq. (|12p picks up significant contributions from the saddle point away 
from the Fermi surface where -^{kx — k 2 t/2m)\k=k* — 0, which implies k* = mx/t. 
In the long-time limit VFt/\x\ — > oo, the saddle point moves to the bottom of the 
band, k* jkp ~ > 0. The corresponding contribution to the hole Green's function 
oscillates with frequency kpjlm and, since the dispersion is parabolic about k = 0, 
it decays as l/yt. The decay is slower than that of the low energy contributions, 
which due to the linear dispersion about ±kp decay as 1/t. This dominant role of the 
saddle point contribution is peculiar to one dimension. In general, in d dimensions, 
the saddle point contribution with k sa would decay as t~ d / 2 , while the Fermi 
surface contribution always decays as 1/t in the noninteracting caseJB 




distance 

Fig. 2. (a) Real part of hole Green's function for free fermions with parabolic dispersion as a 
function of distance x (in units of hp 1 ) and time t (in units of m/fcS.). The dashed line indicates 
the light cone \x\ = vpt. (b) The long time behaviour inside the light cone is governed by high 
energy states at the bottom of the fermion band. The effective model for the interacting case keeps 
states with k fa and A; fa ±kp. 

In the case of parabolic dispersion there is no saddle point above the Fermi level 
for vpt/\x\ — > oo, thus the long time behaviour of the particle Green's function 
G p (x,t) = (*(a;,0* t (0,0)) is well described by LL theory© On the other hand, for 



Turning on weak interactions in d > 2, power law decay is replaced by exponential decay due to 
the finite lifetime of quasiparticles in Fermi liquids. 

c The problem would be the short time limit, due to the lack of a natural high energy cutoff for 
particle states. 
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a tight-binding model with nearest neighbour hopping J the dispersion relation is 
e(k) = — 2 J cos k (setting the lattice parameter to 1). In this case the single-particle 
spectrum is bounded from above as well as from below and there is another saddle 
point at the top of the fermion band. Thus the particle Green's function for a lattice 
model also oscillates in time inside the light cone. 

Numerical methods such as the time-dependent density matrix renormalization 
group (tDMRG)p2 and the time-evolving block-decimation (TEBD) algorithm^ 
provide direct information about the real time evolution of one-dimensional sys- 
tems. These methods reveal that interacting models also exhibit high frequency 
oscillations in equilibrium correlation functions inside the light cone vt > \x\, with 
renormalized velocity i)P3 Although there are no exact analytical results for time- 
dependent Green's functions of interacting models, not even for integrable ones, 
the exact solution for the noninteracting case suggests an approximation to de- 
scribe the long time decay in the interacting case. Besides the low energy chiral 
components, it is important to consider modes with parabolic dispersion at the 
bottom of a renormalized fermion band. This can be achieved by pushing the x-ray 
edge methods for nonlinear Luttinger liquids^ beyond the low energy regimeP^ 
Starting from the non-interacting band, one expands the fermion held in the form 
- e ikFX ip R (x) + e- ikFX ip L (x) + cft(x), where the high energy held cft(x) is 
defined from hole states with momentum k m (see Fig. 2(b)). Using this mode 
expansion for the kinetic energy and interaction terms in the Hamiltonian leads to 
the familiar model of a mobile impurity in a Luttinger liquid The impurity, 
defined in momentum space in this case, is a single deep hole at the bottom of the 
band. The long time limit of the hole Green's function in the interacting case is 
controlled by decay of the deep hole due not only to the parabolic dispersion, but 
also to scattering by low energy particle-hole pairs. The coupling between the deep 
hole and the low energy modes can be treated exactly within the effective impurity 
model using a canonical transformation that shifts the bosonic fields.^ The result is 
that in addition to the (sub-leading) CFT terms in Eq. ([5]) the hole Green's function 
for vt^> \x\ » kp 1 has a time-oscillating term 

-iWt+iMx 2 /2t 

G h (x,t) = -. (13) 

Here W and M are the energy and effective mass of the deep hole, respectively. The 
exponent for large t in Gh{vt/\x\ — > oo) ~ e~ tWt /t^ +v differs from the noninteract- 
ing result by an orthogonality catastrophe correction 



1 

2K 



(14) 



where 5 is interpreted as the phase shift of the Fermi surface states due to the 
creation of the deep hole. For a weak short-range density-density interaction Hi nt = 
(1/2) / dxdx'V(x - x')n{x)n(x'), one finds to lowest order^ S a (Vq — Vk F )/vF, 
where is the Fourier transform of the interaction potential V(x). For integrable 
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models, it is possible to extract phase shifts from Bethe ansatz equations and then 
compute the exact exponents for strong interactions.^^ More generally, phase 
shifts can be determined from information about the exact high energy spectrum 
using phenomenological relations 

The long time behaviour of the density-density correlation function \{x,t) = 
(n(x, t)n(0, 0)) also involves high energy modes. For free fermions, \{x,t) factor- 
izes into particle and hole Green's functions. In the case of a parabolic dispersion 
relation, the longest lived particle-hole excitation has total momentum ±/cp and cor- 
responds to a hole at the bottom of band and a particle at either one of the Fermi 
points. With 1/Vt decay for the hole Green's function and 1/t decay for the particle 
Green's function, the density-density correlation function for free fermions oscillates 
with frequency kp /2m and decays as 1 /t 3 / 2 . Again, this should be compared with 
the decay predicted by LL theory. According to Eq. ([9} , there are low energy contri- 
butions with momentum q ss (particle-hole pair around a single Fermi point) and 
q sa ±2kp (particle-hole excitation between the two Fermi points). In the noninter- 
acting case, both contributions decay as 1/t 2 , more rapidly than the high energy 
contribution. Turning on interactions between the fermions, the parameters of the 
dispersion are renormalized and the exponent of the high energy term in x(x, t) is 
modified by x-ray edge type effects. The general decay is of the form 

±ik F x-iWt+iMx 2 /2t 
y(x t) ~ == ==— (15) 

V~t(vt T xY x -V^) 2 (vt ± X )^-V^) 2 ' 

with A, A defined in Eq. Q. In the long time limit, %(f >• \x\/v) ~ e~ iWt /t n with 
exponent!^ 

l + K 1 / S\ 2 

+ 2k{ 1 -^)- (16) 

In this case, the correction to the free fermion exponent is of first order in the 
interaction: rj sa 3/2 — S/tt for V <C vp. The exponent r] is related to the lower 
edge singularity of the dynamical structure factor S(q,cv) for q = hp. The fact that 
rj decreases with an increasing repulsive interaction is manifested in S(kp,uj) as a 
divergence at the lower edge, similar to the effect observed in the quasi-boson peak 
at low energies (see again Fig. 1(c)). 

For a noninteracting lattice model, the longest lived particle-hole excitation is the 
one obtained from the saddle point contribution for both particle and hole Green's 
functions. This is the excitation with total momentum q — tt that has the maximum 
energy allowed for a single particle-hole pair. For simplicity, let us restrict ourselves 
to the particle-hole symmetric case of a half-filled lattice, kp = ir/2. The density- 
density correlation function picks up two factors of e~ lWt / y/i, where W — 2 J is half 
the bandwidth, from the decay of hole and particle with parabolic dispersion. As a 
result, at large times we obtain x(t ^> \x\/vp) ~ e~ l2Wt /t. However, it turns out 
that this contribution is strongly suppressed by repulsive interactions. The reason 
is that the problem of two high energy particles (or a particle and hole) at the 
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threshold of a continuum where there is an inverse square-root divergence in the 
joint density of states — related to the slow l/t decay of xi x it) m rea l space and 
time — is analogous to the exciton problem in one dimension. For arbitrarily weak 
interactions, resonant scattering between the two particles removes the divergence 
at the threshold of the density of states Consequently, the exponent must change 
discontinuously when the fermion interaction is switched on. For the integrable 
model with nearest neighbour repulsion, V(q) — 4JAcosg with < A < 1 (which 
is equivalent to XXZ spin chain with anisotropy parameter A), it is verified^ that 
the l/t decay turns into a 1 /t 2 decay for times t > 1/ J A 2 . In non-integrable models 
the suppression must be even stronger because this contribution is connected with 
the upper threshold of the particle-hole continuum in frequency domain and the 
power law at this threshold is broadened by coupling to the continuum of multiple 
particle-hole pairsP^In real time, this implies an exponential decay for times larger 
than the corresponding decay rate. The conclusion is that also for lattice models 
with repulsive fermion interactions the long time behaviour of x(x, t) is governed 
by the excitation with one single high energy particle (or hole) and described by 
Eq. (15). 

Summarizing this section, at zero temperature correlation functions of critical 
one-dimensional systems with nonlinear dispersion oscillate at large times and decay 
as power laws with non-universal exponents. The exponents depend not only on the 
Luttingcr parameter but also on phase shifts associated with high energy modes 
with parabolic dispersion. Calculating the exact exponents requires information 
about the exact spectrum. There are, however, exceptions where the exponents do 
assume universal values because they are constrained by the high symmetry of the 
modelPS For example, in spin chains with SU (2) symmetry the phase shift S is 
fixed to 6 = 7r/2 and the oscillating term in the time-dependent spin correlation 
function decays as e~ lWt jt independently of details of the interactions (this includes 
nonintegrable models with finite range spin exchange interactions and even the 
Haldane-Shastry modeP^ with 1/r 2 long-range interactions) . On the other hand, at 
the SU(2) symmetric point the Luttinger parameter becomes^ K = 1/2 and the 
staggered (q — 2kp — it) low energy contribution decays as l/t 2K = l/t at large 
times, i.e. with the same exponent as the high energy contribution. Moreover, for 
SU (2) symmetric models one should also expect logarithmic corrections in the long 
time decay due to marginally irrelevant couplings between high energy modes and 
low energy SU(2) currentsP*] 

4. Long time decay at finite temperatures 

Conformal invariance implies that in the TL model the correlation function for a 
field &(x) with conformal dimensions (A + , A_) decays at finite temperatures as^ 

(17) 



nT/v 
sminrT(x/v — t) 



irT/v 
sinh nT(x/v + 1) 
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where T is the temperature. For values of x and vt which are small compared to the 
inverse temperature, |x±ut| <C v/T, one observes the power law decay characteristic 
of zero temperature correlations. For |a;±i;t| ^> v/T, the correlation function decays 
exponentially ~ e -2TrA + T(x/v-t) e -2irA-T(x/v+t) w ^ a thermal correlation length 
t ~ v/T. 

As discussed in the previous section, for models with nonlinear dispersion the 
slowest decaying term in correlation functions involves the propagator of a deep 
hole coupled to low energy modes. In Eqs. fp~3|) and (fT5|) . for example, the long time 
decay is determined by the factor e~ %wt /y/t from the Green's function of the free 
hole with energy W and parabolic dispersion together with the factors from low 
energy chiral fields. The conformal dimensions of the latter are determined after 
the canonical transformation that decouples the deep hole. For T <C W, the main 
effect of thermal fluctuations is to replace the low energy factors in the result from 
the effective impurity model by the corresponding finite temperature expressions 
according to Eq. ()17[) . As a result, the oscillating term in time-dependent correlation 
functions also decays exponentially within a time scale ~ 1/T. For instance, for the 
density-density correlation function one gets 

X (t » T-\ x/v) ~ (18) 

The exponential decay in real time is connected with the thermal broadening of x- 
ray edge singularities in frequency domain. Beyond the expression in Eq. (1181) . one 
should also include the decay of the high energy hole due to three-body scattering 
processes, which lead to relaxation time Th ~ W/T 2 W^ 

However, the long time behaviour of low energy correlators can also be affected 
by irrelevant operators. Since irrelevant operators introduce interactions between 
bosonic modes, the interesting possibility is that inelastic collisions lead to diffu- 
sive behavior in Luttinger liquids. The word diffusion is used here in the sense of 
phenomenological theories^ for many-body systems at high temperatures which 
predict that the autocorrelation for the density of a globally conserved quantity 
decay as l/t d / 2 in d dimensions due to scattering-dominated random walk of the 
excitations. In one dimension, this means a \/\[t decay, which is clearly much slower 
than the exponential decay predicted by scaling form in Eq. (|17p . 

A finite boson lifetime that gives rise to diffusion is indeed found for lattice mod- 
els at lialf-filling SUH3 At half-filling, particle-hole symmetry rules out cubic band 
curvature terms. In this case the leading perturbations to the TL model are quartic 
band curvature operators, which have scaling dimension four, and the non-oscillating 
umklapp term, which has scaling dimension AK . For the integrable model with near- 
est neighbour interactions only, the coupling constants for these perturbations are 
known exactly! 4 ^ While the calculation of the boson self-energy at zero tempera- 
ture is again plagued by on-shell singularities, at finite temperature the calculation 
is well behaved in the regime \w ± vq\ <C T. The result 42 is that band curvature 
operators only contribute to the real part of the self-energy. The finite temperature 
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Fig. 3. Schematic time decay of density autocorrelation function = { n ( x > t)n(x, 0)) for weakly 
interacting lattice model at half filling and low temperature T. For t <C 1/T, x(t) oscillates and 
decays as a power law. For t > 1/T, the high frequency contribution decays exponentially. For 
t S> l/f{T), the low energy contribution becomes diffusive, xifi) ~ ^vT^OA' 



decay rate is due entirely to umklapp scattering SH U — A f dx cos(4v 7r K cf>) . Com- 
puting the self-energy to first order in band curvature and second order in umklapp, 
the long wavelength part of the density-density correlation function (i.e. the boson 
propagator) becomeqj 

Here A(T) w vK/[l + b(T)} and v{T) w v[l + c(T) - fe(T)] 1 / 2 with b(T),c{T) - 
0(T 2 , T SK ~ 4 ) determined by the real part of the boson self-energy. The decay rate 
is given by32HH 

7 (T) « A 2 ^bcos 2 (2^)r 2 (2^)r 2 Q - 2Jf) T 8K ~ 3 . (20) 

The imaginary part of Eq. (fl9)) yields a Lorentzian peak with width j(T) <C T 1 for 
all K > 1/2. Particularly at the value K = 1/2, for which the umklapp operator 
becomes marginal, the decay rate picks up logarithmic corrections and we obtain 
j(T) ~ T/\d 2 {W/T). We stress that this Lorentzian approximation is only valid in 
the regime where j(T) dominates over the T = broadening due to band curvature 
for finite q. Using Eq. (fT9j) to calculate the correlation function in real space and 
time, one finds that inside the light cone t > \x\/v there is an additional contribution 



d This expression assumes that the model is far from an integrable point and there is no ballistic 
channel contribution to x(?,^)P^ 
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to x( x >t) besides the standard CFT terms. For times t ^> l/-f(T) ^> \x\/v, this 
contribution can be calculated analytically and reads^ 

This is precisely the decay expected for classical diffusion in one dimension, 

even though the assumptions of phenomenological theories definitely do not hold 
for Luttinger liquids. As a result, the long time behaviour of the density-density 
correlation function at finite temperatures is governed by low energy modes, since 
the high energy terms die out exponentially for t > 1/T. Notice that the important 
effect here is the correction to scaling due to the irrelevant umklapp operator, not 
to band curvature operators that break Lorentz invariance. The real time decay of 
the autocorrelation function x( x — 0>*) f° r T ^ is illustrated in Fig. 3. 

The mechanism of diffusion in the density correlator is important to explain 
the spin-lattice relaxation rate observed in spin- 1/2 chains. 45 It is also connected 
with the question of ballistic versus diffusive transport in integrable one-dimensional 
systems P^l 



5. Conclusion and outlook 

We have revisited time-dependent correlation functions of critical one-dimensional 
systems in light of recent theoretical advances that have taken us beyond the 
paradigms of LL theory. Although very successful in describing thermodynamic 
properties, LL theory breaks down in the calculation of dynamical properties in the 
presence of band curvature. When treated as perturbations to the Lorentz-invariant 
TL model, formally irrelevant band curvature operators generate divergent correc- 
tions to the boson propagator. However, it is possible to resum the divergences 
in perturbation theory by refermionizing the bosonic excitations and considering 
particle-hole pairs in a band with renormalized nonlinear dispersion. This proce- 
dure yields a line shape for the "quasi-boson" peak in the dynamic density-density 
response which is remarkably different than the quasi-particle peak in the spectral 
function of a Fermi liquid. 

Also due to band curvature effects, LL theory misses the leading terms in the 
asymptotic long time behaviour of correlation functions at zero temperature. In one 
dimension, high energy modes with parabolic dispersion located at the band edges 
give contributions to time-dependent correlation functions which oscillate in time 
and decay as power laws with smaller exponents than the standard contributions 
from low energy modes. 

At finite temperatures, low energy and high energy contributions in general 
decay exponentially. However, for lattice models at half filling inelastic umklapp 
scattering can give rise to diffusive behaviour of the low energy modes. When this 
happens, the density-density correlation function decays in time as a power law with 
universal exponent ~ 1/Vi for times much larger than the relaxation time. 
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It should be made clear that experimental probes which are only sensitive to 
low frequencies, such as measurements of the local density of states close to the 
Fermi level, p(ui) ~ not affected by the high energy contribu- 

tions discussed here. However, the effects in real time evolution could be observed 
in experiments with ultra cold atoms in optical lattices. Coherent equilibrium dy- 
namics can be investigated by preparing the system in the ground state and then 
creating a local perturbation, which is possible with the development of techniques 
to address individual atoms.^Time oscillations similar to the ones we discussed are 
also seen in numerical simulations of non-equilibrium dynamics of one-dimensional 
systems.^ The oscillatory behaviour is beyond the light cone effect predicted by 
conformal field theory methods for quantum quenches,^ but is presumably inter- 
preted in terms of effective band edges for highly excited states. Experiments suggest 
that questions about the decay of correlations in equilibrium and non-equilibrium 
dynamics require a better understanding of the role played by the integrability of 
the model!^ 

The prediction of a diffusive long time tail at finite temperature is consistent 
with nuclear magnetic resonance experiments which probe the dynamics of spin- 
1/2 chainsP^ The decay rate j(T) has been confirmed numerically through the 
decay of the current-current correlation functionp^J but not directly in the density- 
density correlation function in the low temperature regime. It would be interesting 
to investigate the effects of spin diffusion in other correlation functions, such as 
the single-particle Green's function for spin-1/2 fermions. The study of finite tem- 
perature dynamics of Luttinger liquids should benefit from the recent progress in 
numerical methodsP^E^ 
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